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SUMMARY 


Theoretical  probability  distributions  for  two  detector  statistics 
are  derived  and  used  to  develop  performance  comparisons.  "Statistic  A" 
is  a sum  of  terms  representing  cross-correlations  between  a reference 
channel  and  two  auxiliary  channels  independent  of  the  reference  and 
of  each  other  except  for  a common  sinusoidal  signal  component. 
"Statistic  B"  is  a version  of  statistic  A normalized  by  the  envelope 
of  the  reference  channel.  Which  of  these  statistics  yields  the 
better  detector  is  shown  to  depend  upon  the  ratio  of  the  noise 
power  received  on  the  auxiliary  channels  to  that  on  the  reference: 
when  this  ratio  is  less  than  0.25,  the  normalized  statistic  is  to 
be  preferred.  Computational  procedures  are  fully  disclosed. 
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DISTRIBUTIONS  FOR  TWO  CROSS-CORRELATION 
DETECTOR  STATISTICS 


INTRODUCTION 


A three-channel  conununications  system  is  postulated  to  have  the 
following  elements;  a reference  channel  and  two  auxiliary  channels, 
all  occupying  the  same  bandwidth.  Let  the  reference  channel  output 
be  designated  x (t)  and  the  outputs  of  the  auxiliary  channels,  x (t) 
and  X (t) . For®no  signal,  each  channel  output  is  assumed  to  be  ^ 
independent,  zero-mean  Gaussian  with  the  variances 


= N (reference) 

® (1) 
= aN  (auxiliaries) 

1 2 


The  signal  power  received  on  these  channels  is  taken  to  be 


with 


S = (reference) 

0 

S = 'sk^S^  (auxiliary  #1) 

1 1 

S = ’sk^S^  (au:<iliary  #2) 

2 2 


k^  + k^  = 1 

I 2 


(la) 


(lb) 


In  this  study,  the  probability  distributions  are  derived  for 
two  detector  statistics  obtained  from  samples  of  the  channel  outputs, 
and  some  calculations  of  their  performance  are  presented.  Given 
the  glossary  of  notation  shown  in  Table  1 and  the  detector  model  of 
Figure  l,the  two  statistics  to  be  studied  may  be  written 

1 1/2 
z.  = •=-  Z [Z^  COS^  (())  -(}>)+  Z^  cos^  ((|)  -()>  )]  (2) 

A 2 0 1 1 0 2 2 0 

and 

z = i [z*  cos^  (<!)  -4>  ) + cos^  (<1)  -<(>  )1  (3) 

B 2 1 1 0 2 2 0 
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Table  1 

Glossary  of  Notation 


Subscripts  (unle 

Z 

z 

a 

b 

<P 

k , k 
1 2 

K 

n 

s 

N 

a 

T 

M=WT 

h^=sV2N 

ro 

F (;;) 

1 1 

Q(x) 

Q(x  [v) 

He^ (x) 

4>  (x) 


ss  otherwise  noted):  Oireference;  Itauxiliary  1; 

2: auxiliary  2 

noise  envelope 
detector  statistic 
noise  in-phase  component 
noise  quadrature  component 
phase  angle 

auxiliary  channel  weighting  factors 
signal  in-phase  component 
signal  quadrature  component 
signal  envelope 

noise  variance  (reference  channel) 
auxiliary/reference  noise  power  ratio 
threshold  value 

number  of  post-integrations  (=bandwidth-time  product) 
signal-to-noise  ratio 
gamma  function 

Bessel  function,  third  kind,  integer  order  n 
Confluent  hypergeometric  function 
Gaussian  probability  integral 
Chi-squared  probability  integral 
Herm.ite  polynomial 

Gaussian  probability  density  function 
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MODEL 


REF.  W 


BPF 


AUX.  # 1 


AUX.  # 2 


BPF 


BPF 


SAMPLING  AND 
ARITHMETIC 
UNIT  (MEMORY 
LESS) 


M=WT 


z. 

1 ^ 

M 

Z . 

THRFfsHni  n 

l/M^  2. 

DECISION 

i=1 

z>r 

SIGNAL 

NO 

SIGNAL 

Z<T 


BANDPASS  DETECTOR  POST-INTEGRATION 

FILTERS  ALGORITHM 


INPUTS;  GAUSSIAN  NOISE  PLUS  DETERMINISTIC  SIGNAL 
Xq(x)=Zq(x)  COS(oJ^t-0o)-i-S(t)  COS(co^t-0) 

= [a(,(t)+^(t)]  COS  a;^t+[bQ(t)+T)(t)]  SINcj^t 
X^  (t)=Z^  (t)  COS(OJ^t-  <t>^  l+k,  S(t)  COS(OJ^t-0) 

= [a^  (t)-t-  k.,  ^(t)  ] COSco^t+[b^(t)+k.,  i?(t)  ]siNoj^t 

X2(t)=Z2(t)  COS(oJ^t-02)+k2  S(t)  COS(CO^t-0) 

= [a2(t)+k2  ^(t)]  COS  OJ^t-^[b2(t)+k2  17(1)  ] SIN  W t 
WITH  k.,2  + 1(2^  = 1 

GAUSSIAN  COMPONENTS  ARE  ZERO-MEAN,  WITH 

e|x||.n  I 

> NO  SIGNAL 
Ejx2|=Ejx2|=aN) 


DETECTOR  STATISTICS  (WT=1=M) 

Z^=-\l2yJz\  Z2  COS2(0,-0q)+z2  z|COS2(02-0o) 

Zb=1/2  [z2  COS2(0^-0jj)-hZ2cOS2(02-0o)] 

FIGURE  1 DETECTOR  MODEL 


I 

J 
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Statistic  z can  be  interpreted  as  a normalized  version  of  statistic 
and  botn  represent  combinations  of  correlations  between  reference 
and  auxiliary  channel  outputs.  Successive  calculations  of  z^  and  z^ 
are  considered  to  be  independent. 

The  joint  probability  density  function  (pdf)  of  the  six  in-phase 
and  quadrature  channel  output  sample  components  is  given  by 


p (a  ,b  ,a  ,b  ,a  ,b  ) 
10  0 112  2 


(4) 


= [a^(2iTN)M~^  exp  [a(a  - + a (b  -ri)^ 


2aN 

+ (a  -5k  + (b  - nk  )^  + (a  -5  k )"  + (b  -nk  ) M }, 

1 1 1 1 2 2 2 2 

where 

5 = S cos  <[) 

= S sin  <() 

are  the  signal  components. 

DISTRIBUTION  OF  STATISTIC  A FGR  WT=1 

In  reference  to  the  pdf  shown  in  Equation  (4) , we  define  the 
transformation  of  variables: 


(5) 


a 

0 

z 

0 

cos 

y 

N 

o 

O 

b 

0 

= 

z 

0 

sin 

Y 

0 < Y £ 

a 

1 

X 

1 

cos 

y-  y^  sin  y 

— co<x  <00 
1 

b 

1 

= 

X 

1 

siny 

+ 

y cos  Y 

— oo<y  <00 

1 

a 

2 

= 

X 

2 

cos 

Y 

- y sin  Y 

2 

— oo<x  <00 
2 

b 

2 

ST 

sin 

Y 

+ y cos  Y 

2 

— oo<y  <00 
2 

(6) 


The  Jacobean  of  the  transformation  is  z ^ so  that  the  joint  pdf  of 
the  new  variables  is  ” 

p ( , Yf  X , y , X , y ) = 

2 0 112  2 


2 p (2  cosy,  z sinYfX  cosy-y  siny,x  siny  + y cosy, 
010  0 1 1 1 1 
x^ cosy-y^ siny , x^sinY+  y^cosy) . 


(7) 
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After  integrating  out  y and  y we  obtain 

1 2 


p (Z  , Yf  X , X ) = 


3 0 


1 2 


^ — exp  {- la(Z  cosy- 5) 

ci(2ttN)2  2aN  “ 


(8) 


+ ot(Z  siny-ri)^  + {x  -S  k cos(y~4>)}^  + {x  -S  k cos  (y~(}))  } 

0 11  2 2 


We  now  define  another  transformation: 
Z = V /2z 


Y = ij;  + <() 

X = /2z  cos  c/v 


X 2 = sin  t,/v 


with 


z ^ 0 
V ^ 0 

0 <_  £ 271 

0 < C < 277 


(9) 


The  Jacobean  of  this  transformation  is  2 , so  that  the  pdf 

becomes 


P (z,  V)  = ; 

**  av  (ttN)  ^ 


exp  {-  [a  (2zv^-2/^Sv  cos  iji  + S^) 

2aN 


+ 2z/v^  - 2/Jz  S cos(?-e)  cos  \p/v  + cos^i^;]  } 


(10) 


where 


0 = tan"*  (k  /k  ) (10a) 

2 1 

Using  h^  = S^/2N,  expanding  the  second  and  fifth  exponential  terms 
in  series,  and  eliminating  v by  integration,  we  find 


P (z,  i|),  C)  = exp  {-h^  (1  + cos^  \p/a)  } x 

5 a(7rN)^ 
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^ E E 

m=0  n=0 


1 


mini 


m 

[^h  cos  ^ /yj  ji^ 


cos  (C-6)  cos  ip 


ll" 

nJ 


X 


00 

/-  m-n-l  r 2 

dv  V exp  ^ 

o 


y v"^ } 


exp  {-h^  (l+cos^ijj/a)  } 

a (ttN)  ^ 


00  00 


in=0  n=o 


(11) 


X 


cos  (C-6 ) 


COSlJj 


/I]"  g<.-n,/. 


^(m-n)/2 


f 


where  we  have  used  y “ and  y = — in  integration  formula  3,471.9 
of  reference  1.  Aft4r  writing^ thi^factor  exp{-h^cos^il//a}  in  series 
form  and  using  the  integration  formula  [reference  1,  #3,661.2] 


k 

(a  sin  u + b cos  u) 


0,  k=2r+l 


(12) 


we  successively  eliminate  ? and  \p  by 

(after  recognizing  a series  for  F ) 

1 1 


inegration  to  obtain  the  pdf 


p (z) 
6 


4z  -h^ 
e 

aN^ 


Z Z 
m=:0  n=0 


(h^z/N)"^'^’^  ^-(m+3n)/2  p 

ml  (m+n) ! (n!)^  T (m+H) 


(13) 


X (2z/N  /a)  (m+n+ h ; m+n+l;-h^/a)  . 

The  argument  of  the  pdf  is 

z = Z »^^+x^ 

2 0 1 2 

'l.  S.  Gradshteyn  and  I.  W.  Ryzhik,  Table  of  Integrals,  Series,  and 
Products  (4th  Ed.),  Academic  Press,  New  York,  1965. 

9 


NSWC/WOL  TR  78-37 


= 4 z /o 


^ « r (a  cosy+b  sin  ) ■^  + (a  cosY+b  sinv) 
^01  1 2 2 


= h Z /z^cos^  {Y“4>  ) + Z 


(14) 


^ fc*  r u \ y - sp  / 1 "COS  ^ Ik 

^01  12  2 A 


where  Y - ^nd  M=WT=1.  As  a check,  we  note  that  for  no  signal  (h^=0)  , 
0 


00  00 
J"  dz  (z;h^=0)  = J'  dx  x (x)  = 1. 


(15) 


MOMENTS 


Writing  (13)  in  shortened  form. 


00  00 


P (z)  = 


aN^ 


£ / 1-2  x->\  m+n+1  „ 

f(m,n;a,h  ,N)  z K, 


m=0  n=0 


(n/s) 


' (16) 


and  using  the  integral  [reference  1,  #6.561.16] 

f dx  K (ax) 

J m-n 


m+n+y 


r(m  + I + l)r(n  + I + 1)  , 


(17) 


we  find  that 


00  00 


E{z^}  = (N/a)^  ^ I f (m,n;a,h^  ,N)  (N/a')"'''’^r  (m+^1)  T (n+^1)  . 

m=0  n=0 


For  h^  = 0 , the  moments  are,  explicitly, 
E{z^|h*=0}  = (N/a)"  [r(  1)]^. 


2 ' '“2^ 
(18) 

(19) 
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Thus  for  noise  only, 

( E{z}  = Nir/a74  = .7854N/a 

* a = N/a  = .6190N/cr. 

z 

For  small  h*  we  have,  approximately, 

E{z^}  = (Nv^)^  [ r(  I + 1)  1^1  + ^ (1  + 

PROBABILITY  INTEGRAL 


(20) 


(21) 


Given  the  pdf  (13)  , the  corresponding  (complementary)  probeibility 
integral  is 


Qa(x)  = 


o 

= p^{z  > t}  = y* 


dz  p (z) 
6 


-h 


00  00 


= e I I iiSiili W_L<I>i2±«_  F .hV«) 

m=0  n=0  m!  (n! ) ^ (m+n)  ! r (n+Js) 


X y*  dx  X 


m+n+1 


„(x)  , 

m-n 


(22) 


where 


= 2T/N/a  . (23) 

To  solve  the  integral,  we  use  the  fact  that  K = K_^  and  make  the 
change  of  variable 


X = T /w^  + 1 
1 


to  get 


^ m+n+1  y dw  w (/w^+ 

1 J 


TTTxm+n 


1)“'"“  K „(T  /w^) 
n-m  1 


wK  (t  /w*+  1) 
n-m  1 


(/w*  + 1) 


n-m 


(24) 
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= T 


m+n+1 


I 0 / 


dw 


^2k+lK  /w^+  l) 

n— in  1 


(/vr  + 1) 


n-m 


(25) 


= T 


m+n+1 


where  integral  6.659.3  of  reference  (1)  was  used.  Substituting 

this  result  in  (22),  using  Rummer's  transformation  on  the  hypergeometric 

function  , and  manipulating  indices  gives  us 


Q-  (t)  = t e 
A 1 


-h^ (1+1/a) 


" ( ^h  ^ T ) 


I '~Cm'l ) ^ r (m+?5)  ^F^(is;m+l;h‘ 


/a) 


< I H 

n=0  V"/ 


(2/aT^ ) 


m=0 
n n 


-tr  ^ 
n+^s)  k=«0 


(T  /2) 


(26) 


k! 


^m-n-k+1 1 ^ . 


For  noise  only,  is  the  probability  of  false  alarm  and  is  simply 

Q (Tlh2=0)  = T K (T  ) = -^  Ki  1-^  .-7, 

' N/~  \N/ay*  ^27) 

DISTRIBUTION  OF  STATISTIC  B FOR  WT  = 1 

Beginning  with  p (z  , y#  x , x ) as  written  in  equation  (8),  we 

3 ® 12 

make  the  following  transformation  of  variables: 

0 < X < " 

with  0 <_  6 ^211 

Y=\{)  + <j)  0<i)j<2tt 


X = /2X  cos6 
1 

X = /2X  sin6 
2 


(28) 


The  Jacobean  of  the  transformation  is  1,  giving 

p^(Z^,  X,  5,  exp{-  -^[a  (Z*-2Z^  ScosiJ.+s' ) 


+2X-2S/2Xcosijicos  (6-0)  +S^cos*ij)]  } 


12 


(29) 
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Eliminating  by  integration,  we  get 

-h^  -1 

V (X,  S , ip,)  = — - e"  exp{  :~[2X-2S/2rcosi);cos  (6-6) 

■^8  aN(2TT)^  2aN 

+ S^cos^il^]} 

X ! r(  I*  1),  using  = |i. 


Integrating  over  6 results  in 


4 > 


(2hcosT;')  _ / S . 1 \ 
ml  1 V 2 


Expanding  the  Bessel  function  in  its  series  form  and  using  (12) 
to  integrate  over  ij;,  we  obtain  2n 

- 1 ^ r 7 1^) 

“ aN  ■ aN^  [ r(m+J5)  ^ , , 

m=0  n=0  nin! 


(-h^/g)  r (m+n+k+>s) 
k ! (m+n+k) ! 


m=0  n=0 


By  using  Kximmer's  transformation 


j^Fj^(a;  b;  -c)  = e~  (b-a ;b ; c)  , 

on  (32) , an  alternative  expression  for  (32)  is  found  to  be 
1 


1 ^ f v,2mxi/  ^ V 7 h^“‘(X/Na^) 

aN  (1+1/a)  - ml  (nl)^ 


nmtis) 


^Fj^(i5;  mfl;  h*/a) , 
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The  argument  of  this  pdf  is 

X = 35  (x^  + x^) 

1 2 

= islZ^cos^  (Y-<|)  ) + Z^cos^(Y-<j>  )]  = z (35) 

1 1 2 2d 

where  y - <|)  and  WT  = 1, 
o 

When  h*  = 0,  (32)  and  (34)  reduce  to 

^io(X;h*=0)  = ~ exp{-X/aN} , (36) 

the  exponential  distribution  with  parameter  aN.  Therefore,  to  a 
factor,  (for  no  signal)  is  distributed  as  a chi-squared  variable 

with  two  degrees  of  freedom; 

2Z 

^ X"  (2)  for  h*=  0.  (37) 


MOMENTS 


Using  (34) , we  find  that 


E{X^}  = (Nct)^  Vl 


“ m 

^ I I im+h)  (y+Pn  i .^2 npl 

m=0  n=0  m!  (n!  ) '^r  (m-n+ij)  (35,m+l,h  /a).  (38) 

For  small  h^ , this  expression  reduces  to 

E{x’'^}  = (Na)’"^  y!(l+  hV2a);  (39) 

thus  for  h*  = 0,  both  the  mean  and  standard  deviation  equal 
Na — a result  which  is  predictable  in  view  of  (37) . 

PROBABILITY  INTEGRAL 

Again  using  (34),  the  (complementary)  probability  integral 
which  corresponds  is  found  from 

Qb(t)  = Pj.{A  ^ xl 


^ g-hM  1+1/a) 


<*>  m 

I I 


h2"'a"^r  (m+>2) 


m=0  n=0  ml  (n! ) * r (m-n+>5) 

00 

J"  du  u”  e = x/aN, 


j^Fj^(i5;m+ljh*/a) 


(40) 
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in  which  the  integral  is  the  incomplete  gamma  function  r(n+l;T^), 
and  is  related  to  the  chi-squared  probability  integral  Q(x*|v); 


_ n . 

r(n+l;T  ) = n!Q(2T  | 2n+2)  = n!  e"  i I rVk! 

^ k=0  ^ 

Putting  the  last  expression  in  (41)  into  (40),  we  have 

^ ~ m n,2m-n,,„.k 

. exp{  - 3^  -h'd+l/c))  I I I ° 

m=0  n=0  k=0 


m!  n!  k! 


(41) 


RECEIVER  OPERATING  CHARACTERISTICS  FOR  OT  = 1 


(42) 


In  order  to  compare  the  properties  of  statistics  A and  B, 
we  shall  first  obtain  receiver  operating  characterstics  from  the 
expressions  already  derived,  that  is,  for  WT  = 1.  in  the  next 
section,  an  approximation  technique  will  be  used  to  extend  the 
comparison  to  arbitrary  WT. 


The  relationships  known  as  receiver  operating  characteristics 
may  be  expressed  in  the  present  application  by  (h’’'  = SNR) 

Pjj  = f(h^•  Pp^,  a,  WT)  (43) 

in  which  the  decision  model  depicted  in  Figure  (1)  is  assumed; 


Pd  - pr(2  > T^)  . i = A,B 

where  is  the  false  alarm  threshold  determined  from 

a 

^FA  ^ ®i^^a'  “ Pr{z  > x^;  h^=  O}. 

For  statistic  A,  from  (27)  we  have 


'’fa  = \a  ''l(^a)'  \a  = 2x^/N/a. 

Using  Table  9.8  of  reference  (2),  the  following  values  can  be 
calculated; 


'’fa 

X 

1 a 

0097=10"* 

5.8 

00103=10"® 

8.2 

000105=10"'* 

10.6 

and  I.  A.  Stegun, 

eds . , 

(Statistic  A) 
WT=1 
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For  statistic  B,  from  (42)  it  is  evident  that 

= exp{-T/aN},  (48) 

from  which  we  may  calculate 


P 


FA 


T /Na 

Si 


0.01 

0.001 

0.0001 


4.60517 

6.90776  (Statistic  B)  (49) 

9.21034  OT  = 1 


Calculation  of  the  probability  integrals  to  find  probabilities 
of  detection  is  more  involved,  but  relatively  straightforward. 

Using  the  computational  approaches  and  programs  described  in  the 
appendices,  (equation  (26))  and  (equation  (42))  were  obtained 

for  several  values  of  the  auxiliary-to- reference  power  ratio  a,  as 
shown -'in  Figure  (2). 

As  expected  an  obvious  feature  of  the  information  displayed 
in  Figure  (2)  is  that,  for  fixed  reference  power  (N),  performance 
is  improved  (smaller  SNR  required)  as  the  noise  power  on  the 
auxiliary  channels  is  decreased  (a  decreased) . What  is  interesting 
in  this  figure  is  what  it  reveals  about  the  relative  performance 
of  the  two  detector  statistics.  It  is  evident  that,  because  of 
different  sensitivities  to  a,  statistic  A is  better  for  a > 0.25 
but  statistic  B is  better  for  a ^ 0.25.  The  interpretation  seems 
to  be  that  the  reference's  amplitude  information  (Z^)  becomes  more 

important  as  the  auxiliary  channels  become  more  noisy. 

PERFORMANCE  STUDY  FOR  W^r  > 1 

To  illustrate  the  extension  of  the  theoretical  results 
we  have  obtained  for  WT  = 1 to  cases  in  which  V7T  is  arbitrary, 
we  examine  the  crossover  effect  of  a on  the  relative  performance 
of  the  two  statistics  for  consistency  as  WT  increases.  Specifically, 
we  compare  values  of  minimxam  detectabxe  signal  (MDS)  for  P = 0.5  and 
Pp^  = 0.01,  expressible  as 

KDS  = hMa,  WT;  P_  = .5,  P„  = .01}  (50) 

D ' FA 

Using  the  moments  for  WT  shown  in  (18)  and  (38)  as  inputs 
to  the  Cornish-Fisher  expansion  described  in  Appendix  B,  we  find 
MDS  by  varying  the  SNR  until 

T(h*;  Q = .5,  a,  WT)  = T(h*=0;  Q = .01,  a,  WT)  (51) 

where  t is  the  inverse  function  of  Q(t) , the  complementary  probability 
integral.  The  results  of  this  procedure  are  shown  in  Figure  (3) , in 
which  we  see  that  a = .25  continues  to  be  a crossover  value  for 
choosing  the  better  statistic. 
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FIGURE  2 RECEIVER  OPERATING  CHARACTERISTICS  (WT=1),  NOISE  POWER  RATIO  (G)  VARIED 
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The  behavior  of  the  MDS  as  shown  in  Figure  (3)  is  so  consistent 
that  it  gives  us  confidence  in  further  extending  the  results  of 
Figure  (2)  with  such  approximate  generalizations  as 


MDS(WT)  = MDS(WT=1)  - .5  - 5 log^^^  (WT)  , 

where  MDS  is  specified  in  dB. 


(52) 


Taking  the  data  of  Figure  (3)  and  making  VIT  the  family 
parameter  results  in  Figure  (4) . In  this  presentation,  the  crossover 
behavior  we  have  been  noting  is  displayed  directly.  The  curvature 
of  the  results  for  statistic  A discourages  us  from  making  any  bold 
statements  about  the  general  dependence  upon  a.  However,  if  we 
make  use  of  the  central  limit  theorem  to  say,  as  WT  -*■  ». 

z ~ N(m2,,  ) (53) 

where  m and  are  the  mean  and  variance  for  WT  = 1,  then  asymptotically 
(51)  becomes 


In 


m (h^)  = m (h^=  0) 
z z 


addition,  since  h^<<l. 


r 9m 

r- 

o 

II 

£ 

z 

9h* 

- 

+ X a (h^=0)//OT. 
we  may  write 
x^Oi  =0)//WT. 


(54) 

(55) 


Finally,  if  as  in  reference  3 we  define 


AMDS  = hV(x 

9m 

= o^(h  =0)/[  — ^ 
^ 9h2 


(56) 


we  can  compare  asymptotic  behavior  of  the  two  statistics  on  a common 
basis . 


From  (20),  (21),  and  (39),  we  get 
AMDS^  = 1.576/(l+l/2a) 

AMDS„  = 2a. 

ID 


(57) 


Figure  5,  a plot  of  (57),  confirms  what  we  have  already  concluded 
from  the  previous  results,  namely,  that  statistic  A provides  better 
detection  performance  for  a > .25,  approximately.  In  addition, 
for  a>>l,  statistic  A's  MDS  relative  to  the  baseline  (which  happens 
to  be  that  of  a square  law  statistic)  is  no  worse  than  2 dB.  We 
see  also,  that  although  B is  the  better  detector  for  a < .25,  A is  at 

^C.  N.  Pryor,  "Calculation  of  the  Minimum  Detecteible  Signal  for  Practi- 
cal Spectrum  Analyzers,"  NOLTR  71-92,  2 Aug  1971. 
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NOISE  POWER  RATIO  (a),  LOG  SCALE 


FIGURE  4 MINIMUM  DETECTABLE  SIGNAL  VS  NOISE  POWER  RATIO, WT  VARIED 
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RELATIVE  MDS  IN  dB 
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2^  2^  1 4 16 


NOISE  POWER  RATIO  (tt),  LOG  SCALE 
FIGURE  5 RELATIVE  MDS  VS  NOISE  POWER  RATIO 
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most  worse  by  2 dB  here,  too;  whereas  when  B is  the  worse  detector, 
its  performance  degrades  indefinitely  as  a increases. 

CONCLUSION 

Further  comparisons  can  be  made  with  the  theoretical  results 
and  computational  procedures  documented  in  this  report. 
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APPENDIX  A 


OUTLINE  OF  COMPUTATIONS 

The  general  form  of  the  computations  which  are 
to  be  discussed  in  themselves  is  the  infinite  series 

f (x)  = f (x)  exp{-h^  -hVct}  S (x) 

1 2 2 

where,  given  the  values  of  the  arguments, 

“ m n 

S = J A(m)G(m)  ^ B(m,  n)  I C(k). 

^ m=0  n=0  k=0 

Truncating  the  series  at  m = and  expanding  on  (A- 
the  computational  approach  we  write 

S = A(0)G(0)B(0,0)C(0) 

2 

m^jo 

+ I A(m)G(m)  {b  (m,  0 ) C (m,  0 , 0 ) 
m=l 

m n 

+ B (m,n)  [C  (m,n,0)  + C(m,n,k)]} 

n=l  k=l 

Common  to  each  of  the  specific  forms  which  are 
below  is  the  function 

G(m)  H rn+1;  h^/a) 

2ma 

“ (2m-l)h^  [(m-l+h^/a)  G (m-1)  - (m-1)  G (m-2)  ] 

G(0)  and  G(l)  are  computed  directly  from 


special  enough 
(A-1) 

(A-2) 

2)  to  show 

(A-3) 

treated  separately 
(A-4) 

, m > 2 . 


A-l 
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F,  (-’5;  1;  X)  = [ ^ 

^ r=0  (r:; 


(h) 


roo 

=•  I D (r) 
r=0 


using  (r)  = 


(rl) 


• ^ 2 


ih) 


v/ith 


D„(0) 


= X (r-Jj)  (r-1) /rS  r ^ 1, 


= 1, 


Similarly,  2;  x)  is  computed  using 

<'=>r 

= FT  TFTTTT 


= x(r  - »5)D^(r-l)/r  (r+1), 
with  D^(0)  = 1. 


(A-5) 


(A-6) 


(A-7) 


Satisfactory  results  were  accomplished  when  truncation  of  the  series 
was  such  that 

D-(r^)  < 10"' ^ (A-8) 

0 “> 

The  form  (A-2) , besides  being  more  convenient  than  a triple 
infinite  summation,  was  chosen  to  imitate  the  power  series 
m“ 

S = ^ (hM"*  E(m)  , (A-9) 

m=0 

with  the  idea  being  tbat  h^  is  most  often  the  major  parameter 
(sometimes  also  being  large)  and  truncctticn  such  that 

(hM"^E(m  )<.0001  S (A-10) 

00  2 

will  insure  convergence  in  h* . For  the  most  part  this  has  worked 
v;ell,  although  for  a>2  and  a<J5  in  some  cases  a more  sophisticated 
approach  (not  assuming  so  uniform  convergence)  might  be  warranted, 
judging  from  the  behavior  of  Q^(t)  and  ©^(t),  for  example,  when 

their  values  are  0.9  and  higher. 
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Q^(t) — Equation  (26) 

f = T 
2 i 

A(m)  H (h^T  /2)™/(m:)* 

' I 

= (h^T  /2)  A(m-l)/m*,  m ^ 1 
A(0)  = 1 

^ (n)  ri^ 

(m-n+1)  (m-n+*i)  B(m,  n-l)/n,  n ^ 1 

B(m,0)  = B(0,0)  = 1 

C(m,n,k)  = C (k)K  ) 

1 m— n— K+i  1 

(k)  = (Ty2)’'/k; 

= {xyi)  (k-l)/k,  k > 1 

C (0)  = 1 

1 

= 2rK^(T^)/T^+K^^(Tj,  r > 2 

T, 

e K (t  ) and  e K (t  ) tabulated  in  reference  (2) . 

0 1 11 


(A-11) 


(A-12) 


(A-13) 


(A-14) 


(A-15) 
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Qg(T) — Equation  (42) 

f = exp(-T/aN)  (A- 16) 

2 

A(m)  E (h^)"'/m! 

= (h^)  A(m-l)/m,  m ^ 1 (A-17) 

A(0)  = 1 


B(m,n)  H r (m+Jj) /a'^nl  r (m-n+^s) 

= (m-n+*5)B  (in,n-l)/an,  n ^ 1 (A-18) 

B(m,0)  = B{0,0)  = 1 

C(k)  = (T/aN)Vk: 

= (T/aN)  C(k-l)/k,  k > 1 (A- 19) 

C{0)  = 1 

y 

E{2^} — Equation  (18) 

f i (N/a)^  (A-20) 

2 

A(m)  = (hM’>(m.')  * 

= h*A(m-l)/iti^  , n'  1 1 (A-21) 

A(0)  = 1 


B(m,n)  E 


B(m,0)  s: 
B(0,0)  = 
C (m,n,k) 


r (m+ij)  r (m+1)  r (m-n+u/2+1)  r (n+y/2+1) 
r (m-n+*i)  r (m-n+1)  (nl ) * a” 

(m-n+>5)  (m-n+1)  (n+y/2)  B (m,n-l)  / (in-n+p/2+1)  an*  , 

(m+y/2)B(m-l,0) , m ^ 1 

[r(y/2+l)  ]* 

= 0,  k > 1 

1,  k = 0 


n ^ 1 
(A-22) 


(A-23) 
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E{Zg} — Equation  (38) 

f E (Na)'^  r(y+l)  (A-24) 

2 

A(m);  same  as  (A-17) 

B(m,n)  E (y+1)^  r (m+J5)/a"(n: ) (m-n+ij) 

= (y+n)  (m-n+Jj)  B (m,  n-1) /an^  , n ^ 1 (A-25) 

B(m,0)  = B(O.O)  = 1 
C(m,n,k):  same  as  (A-23) 
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APPENDIX  B 


COMPUTATION  OF  INVERSE  PROBABILITY  INTEGRAL 

An  asymptotic  expansion  for  arbitrary  probability  distributions, 
adapted  from  § 26.2.49  of  reference  (2),  is  useful  for  computing 
minimum  detectable  signal  (MDS) . Known  as  the  Cornish-Fisher 
expansion,  it  can  be  written 

T ~ m + aW(x  )//M  (B-1) 

p p 


with 


w(x)  = X + [yd] 
1 2 


+ [y  d +Y^d  ] 

2 3 1 4 


+ [Yd  + YYd  + Y*d  ] 

3 5 126  17 


+ 


[y  d +Y^d  +Y  Y 

^ 8 2 9 1 


d 

3 10 


+yJy 


d 

2 1 1 


+Y'*d 

1 


(B-2) 


+ • • • 

where,  by  definition 

Pr{z  > T } = p.  (B-3) 

The  random  variable  z,  the  detector  decision  variable,  is  assumed 
to  be  the  mean  of  M independent,  identically-distributed  samples 
{z^}  with  moments  and  cumulants  as  follows: 


M 


i=l 


^i 


m = < = e{ z.  } 
1 1 


a = K = Var{ z . } 
2 2 1 


(B-4) 


(B-5) 
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K =E{Z?}-3kic  -k^ 

3 1 2 1 1 


K = EiZ?}-4K  K - 3<*-6k: 

4 1 3 1 2 2 11 


K = E{z®}-5k:  k - IOk  k -10k  k^-15k  k^-10k  k^-k®  (B-6) 

5 1 41  23  31  12  211 


K = EfZ?}-6K  K -15k  k -10k-15k  k^-60k  k k 
6 1 51  24  3 41  123 


-15k®-20k  k^-45k^k*-15k  k‘*-k® 

2 3 1 1 2 2 1 1 


Y = K / (Mk  ^ 

1 3 2 

Y = K /Mk^ 

2 4 2 

Y = K / (B-7) 

3 5 2 


Y = K /M^ic^ 

4 6 2 

The  number  and  the  coefficients  {d.}  are  related  to  the  Gaussian 
P 1 

distribution  and  are  given  in  Table  (2)  for  several  values  of  p, 
where  columns  2-4  were  taken  from  page  936  of  Reference  (2) . 

The  brackets  around  the  terms  in  (B-6)  correspond  to  orders 
of  magnitude  with  respect  to  M.  A test  case  run  for  the  non-central 
chi-squared  distribution  and  using  only  the  first  two  bracketed 
terms  yielded  results  8%  and  3%  below  the  true  value  of  MDS  for 
M = 5 and  M = 50 , respectively. 

These  coefficients  can  also  be  used  to  approximate  the 
probability  density  function  via  the  Edgeworth  series f as  shown 
in  reference  (4) . 


’’L.  E.  Miller,  "Computing  R.O.C.  for  Quadratic  Detectors,"  NSWC/WOL 
TR  76-148,  10  Oct  1976. 
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TABLE  2 

COEFFICIENTS  FOR  CORNISH-FISHER  EXPANSION 


p 

.0001 

.001 

.01 

.1 

.5 

.9 

3.71902 

3.09022 

2.32635 

1.28155 

0 

-1.28155 

2.13852 

1.42491 

.73532 

.10706 

-.16667 

.10706 

1.67838 

.84331 

.23379 

-.07249 

0 

.07249 

-2.34115 

-1.21025 

-.37634 

.06106 

0 

-.06106 

.92761 

.30746 

-.00152 

-.03464 

.02500 

-.03464 

-5.17267 

-1.  89355 

-.17621 

. 14644 

-.08333 

.14644 

4.87514 

1.86787 

.25195 

-.11629 

.05247 

-.11629 

.35118 

. 04591 

-.03176 

.00227 

0 

-.00227 

-2.62416 

-.59060 

.07888 

.00776 

0 

-.00776 

-3.48080 

-.70464 

.16058 

.01086 

0 

-.01086 

17.56966 

4.29304 

-.32621 

-.10858 

0 

.10858 

-12.61271 

-3.32708 

.07286 

.09585 

0 

-.09585 
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APPENDIX  C 


COMPUTER  PROGRAMS 


BASIC  programs  for  computing  RCC  and  MDS  are  listed  in  the 
figures  which  follow.  Subroutines  common  to  two  or  more  programs 
are  listed  separately.  Given  the  computational  outlines  of 
Appendix  A,  the  listings  are  nearly  self-explanatory.  Additional 
comments ; 


Program  QAS.  (Figure  6).  Given  the  values  of  a, 

and  computes  Q;^(t)  for  the  values  of  h^  specified  by  the 

user. 

Program  QBS«  (Figure  7) . Given  the  values  of  a and  t , 

cl 

computes  Qg(T)  for  the  values  of  h*  specified  by  the  user. 


Program  MPA.  (Figure  8)  . Given  the  values  of  a,  and  \fT , 

computes  false  alarm  threshold  via  Cornish-Fisher . Then  computes 
detection  thresholds  for  the  values  of  h^  specified  by  the  user,  who 
interpolates  to  get  MDS. 

Program  MDB.  (Figure  9) . Same  as  MDA  but  for  statistic  B. 
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PROGRAM  QAS53A 

1 DIM  G(1I0),K(11() 

2 Pt=4*ATN<I) 

5 REN  PROGRAM  TO  COMPUTE 

7 PRINT  "PROBABILITY  INTEGRAL  FOR  STATISTIC  A (UT=t)" 

If  READ  A1,TI,A 
12  DATA  1,5. 8, .01 

15  PRINT  "GIVEN:  ALPHA, TAU1 ,PFA  ="A1,T1,A 
14  PRINT 

17  PRINT  '■H<DB)”,”PD",”LAST  M" 

21  GOSUB  700 
25  K2=-1 
30  H=1ftK2 
35  X=H/A1 
40  N=0 

45  GOSUB  80f 
250  A2=B2=1 
40  C2=1 

45  S2=G(0)*K(1 ) 

71  M=M+1 

75  K(M+1)=2*M»K(N)/T1*K(N-1) 

81  IF  M<2  THEN  70 
85  GOSUB  840 
90  A2=A2»H»TI/2/Mf2 
110  B3:B2 
105  S1=B2*K(M+1) 

110  FOR  N=1  TO  M 

115  B3=B3»2»(M-N+1 )» (M-N+.5)/N/A1 /T1 

120  C2=1 

125  S0=K(M-NH) 

130  FOR  J=1  TO  N 

135  C2=C2*T1/2/J 

140  IF  (M-N-J+1)<0  THEN  155 

145  30=S»+C2«K(M-N-JH) 

150  GOTO  140 

155  S0=S0*C2»K(N+J-M-1 ) 

140  NEXT  J 
145  S1=SUB3*S0 
170  NEXT  N 
175  L=S2 

181  S2=S2+A2»G(M)»S1 
185  L=(S2-L)/S2 
188  L2=l 

190  IF  L>.00f1  THEN  70 
195  P=S2*T1*EXP<-H-X) 

194  IF  P<.5  THEN  202 

197  IF  L2=1  THEN  202 

198  L2=1 

199  L=I0«L 

200  GOTO  190 

202  PRINT  10*K2,P,N 

203  IF  P<L1  THEN  300 

204  L1=P 

205  IF  P>.99  THEN  300 
210  K2=K2+.1 

215  IF  P.>2«A  THEN  30 
220  K2=K2+.4 
225  GOTO  30 
300  STOP 

700  REH  SUBROUTINE  FOR  BESSEL  FUNCTION 
705  K(0)=.510I258I83*EXP(-5.B) 

710  K(1)=.5524474495*EXP(-5.8) 

715  RETURN 

FIGURE  6 PROGRAM  TO  COMPUTE  PROBABILITY  INTEGRAL  FOR  STATISTIC  A (WT=1) 
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PROGRAM  Q6S53A 

5 HIM  G(109) 

10  PI=4»ATN(1) 

15  REM  PROGRAM  TO  COMPUTE 

2»  PRINT  "PROBABILITT  INTEGRAL  FOR  STATISTIC  B <UT=I)" 

25  READ  A1,T,A 

31  DATA  1,4. 60517, .01 

35  PRINT  "GIVEN:  ALPHA,  TAU,  PPA  "A1,T,A 

40  PRINT 

45  PRINT  "H(DB)'',"PD","LAST  M" 

50  PRINT 
55  T1=T 
60  K2=-2 
65  H=10tK2 
70  X=H/A1 
75  M=0 

80  GOSUB  800 
85  A2=B2=1 
95  S2=G(0) 

101  «=H+1 

105  IF  M<2  THEN  115 

110  GOSUB  860 

115  A2=A2»H/M 

125  B3=B2 

130  S1°B2 

135  FOR  N=1  TO  M 

140  B3=B3»(M-Nt.5)/A1/M 

145  C2=S0=1 

150  FOR  J=1  TO  N 

155  C2=C2*T1/J 

160  S0::S0+C2 

165  NEXT  J 

170  S1=SUB3*S0 

175  NEXT  N 

180  L=S2 

185  S2=S2+A2*G(M)*S1 
190  L=(S2-L)/S2 
195  IF  L>.0001  THEN  100 
200  P=S2*EXP(-T1-H-X) 

205  PRINT  10»K2,P,M 
207IF  P<L1  THEN  300 
208  L1=F 

210  IF  P>.99  THEN  300 

215  K2=K2+.1 

220  IF  P>2»A  THEN  65 

225  K2=K2+.4 

230  GOTO  65 

300  STOP 
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PROGRAM  HDA53A 

5 DIM  G(1l0),V(&),E(6),K(i) 

6 DIM  D(I2) 

II  REM  PROGRAM  TO  COMPUTE 

IS  PRINT  "FA  THRESHOLU  AND  NOS  FOR  STATISTIC  A" 

21  READ  At,A,Mt 
2S  DATA  1,.I1,1I 

31  PRINT  "GIVEN:  ALPHA,  QA,  UT  <"A1,A,MI 
35  PRINT 
41  GOSUB  ill 

45  PRINT  "H(DB)",'‘TA","LAST  M" 

51  PRINT 
55  K3=-.2 
57  K4=-.11 
il  K5=.I1 

65  FOR  K2=K3  TO  K4  STEP  K5 

71  H=1ltK2 

75  X=H/A1 

81  M»=2 

35  GOSUB  811 

90  FOR  1=1  TO  6 

95  G=V(I) 

Ilf  N=l 

115  S2=G(0)*6t2 

III  A2=1 
115  B3=Gt2 
121  N=N+I 

125  IF  N<N9  THEN  135 
131  GOSUB  861 
135  A2=A2»H/Ht2 
Ml  B3=B3»(Mfl/2) 

145  S1=B2=B3 
151  FOR  N=1  TO  M 

155  B2=B2MN-N*.5)*<H-N  + 1)»<N+I/2)/(N-H+I/2+l)/A1/Nt2 
lil  S1=SUB2 
US  NEXT  N 
147  L=S2 

170  S2=S2tA2»G(N)*S1 

175  L=(S2-L)/S2 

180  IF  L>.0II1  THEN  120 

185  E(I)=S2»A1t(I/2)*EXP(-H-X) 

190  M9=M 

210  NEXT  I 

205  GOSUB  1110 

210  GOSUB  1500 

215  PRINT  10«K2,T,H9 

220  NEXT  K2 

300  STOP 

601  REN  DEFINE  GANHA(UI/2)=V(I) 

605  V(1)=SQR(ATN(I)) 

610  U(2)<1 

615  V(3)=V(1)*3/2 

620  U(4)=2 

625  V(5)=V(3)*5/2 

630  V(i)=6 

635  REM  FALSE  ALARM  THRESHOLD 

640  FOR  1=1  TO  4 

645  E(I)=V(I)f2*A1t(I/2) 

650  NEXT  I 
655  GOSUB  1100 
640  GOSUB  1500 

645  PRINT  "FALSE  ALARM  THRESHOLD  =‘T 
470  PRINT 
675  A=.5 
680  RETURN 

FIGURE  8 PROGRAM  TO  COMPUTE  FA  THRESHOLD  AND  MDS  FOR  STATISTIC  A 
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PROGRflH  HIIB534 

5 DIH  G(1»0),V<6),E<6),K(6) 

6 DIM  D(12) 

It  REM  PROGRAM  TO  COMPUTE 

15  PRINT  "FA  THRESHOLD  AND  MDS  FOR  STATISTIC  B" 

21  READ  A1,A,M1 
25  DATA  1,.ei,1B 

31  PRINT  'GIVEN:ALPHA,  PFA,  UT  ="AI,A,Mt 
35  PRINT 
At  GOSUB  600 

45  PRINT  "H(DB)","TA”,'’LAST  M" 

50  PRINT 
55  K3=-.2 
57  K4  = -.n 
60  K5=.01 

65  FOR  K2=K3  TO  K4  STEP  K5 

70  H=10tK2 

75  X=H/A1 

80  M9=2 

85  GOSUB  800 

90  FOR  1=1  TO  6 

100  M=0 

105  S2=G(0) 

110  A2=1 
120  M=M+1 

125  IF  M<M9  THEN  135 

130  GOSUB  860 

135  A2=A2«H/M 

145  S1=B2=1 

150  FOR  N=1  TO  M 

155  B2=B2*(I+N)»(M-N+.5)/A1/Nt2 

160  SI=SUB2 

165  NEXT  N 

167  L=S2 

170  S2=S2*A2^G(M)*S1 

175  L=(S2-L)/S2 

180  IF  L>.0001  THEN  120 

185  E(I)=S2»V(I)*A1tI»EXP(-H-X) 

190  M9=M 

200  NEXT  I 

205  GOSUB  1100 

210  GOSUB  1500 

215  PRINT  10»K2,T,M9 

220  NEXT  K2 

301  STOP 

600  REN  DEFINE  FACTORIAL 

605  V(1)=I 

610  V(2)=2 

615  V(3)=6 

620  9(4)=24 

625  V(5)=120 

630  V(6)=720 

635  REM  FALSE  ALARM  THRESHOLD 

640  FOR  1=1  TO  6 

645  E(I)=V(I)*A1fI 

650  NEXT  I 

655  GOSUB  1100 

660  GOSUB  1500 

665  PRINT  "FALSE  ALARM  THRESHOLD  ="T 
670  PRINT 
675  A=.5 
680  RETURN 

FIGURE  9 PROGRAM  TO  COMPUTE  FA  THRESHOLD  AND  MDS  FOR  STATISTIC  B 
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80*  REH  SUBROUTINE  FOR  IF1 ( ;X) 

8»5  A9=B9=U=U=J9=1 

810  A9=A9*X»(J9-.5)/J992 

815  B9=B9*X*(J9-.5)/J9/<J9+1) 

820  U=U+A9 

825  U=V+B9  I 

831  IF  A9<1E-13  THEN  845 
835  J9=J9+1 
84<  SOTO  810 

845  6(f)=U  \ 

851  6(1)=V 
855  RETURN 

860  D(N)=2*M*<(«-1+X)»G(H-1)-(M-1 )*G(H-2))/(2»M-1)/X 
845  IF  G(N)<G(H-1)  THEN  875 
870  G(H)=1 

875  IF  G(N)>=1  THEN  885 
880  G(N)=1 
885  RETURN 

1108  REH  CONFUTE  CUMULftNTS 
1105  K(1)=E(1) 

1110  K(2)=E(2)-K(1)t2 

1115  K(3)=E(3)-3*K(2)*K<1)-K(1)t3 

1 120  K<41=E(4)-4»K(1)»K(3)-3*K(2)t2-6*K(2)*K(1)f2-K(1)t4 

1 125  K(5)=E(5)-5»K(1 )*K(4)-18»K(2)*K(3)-10»K(3)*K(1 >t2 

1130  K(5)=K(5)-15*K(1)*K(2)t2-10*K(2)*K(1 )T3-K(1)t5 

1135  K(4)=E(6)-4*K(1)*K(5)-15*K(2)*K(4)-10*K(3)t2-15*K(4)»K(1)t2 

1140  K(4)=K(4)-40*K(1)*K(2)*K(3)-15*K(2)t3-20*K(3)»K(1)43 

1145  K(4)=K(4)-45*K(1)t2»K(2)t2-15*K(2)*K(1)94-K(1)t6 

1148  H2=S0R(N1) 

1150  R1=K(3)/H2/K(2)t1.5 
1155  R2=K(4)/N1/K(2)f2 
1160  R3=K(5)/H2t3/K(2)f2.5 
1165  R4  = K(4)/HH2/K(2)t3 
1170  RETURN 
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150«  REM  CORNISH  FISHER  RflUIINE  FOR  IHVERSE  PROfiAOUinr  IHTEG. 

1S0S  REM  AHHISSIBLE  PROS:  .001, .01, .5 

1510  IF  A<.009  THEN  1550 

1515  IF  A>.1  THEN  1580 

1520  0(11=2.32435 

1522  0(2)=. 73532 

1 524  Ii(3)  = . 23379 

1525  D(4)=-. 37434 

1524  11(5)=-. 00152 

1528  11(4)  = -.  17421 

1530  D(7)=. 25195 

1532  0(8)=-. 03174 

1534  0(9)=. 07888 

1535  0(10)=. 14053 
1534  0(11)=-. 32421 
1540  0(12)=. 07284 
1545  GOTO  1405 
1550  0(11=3.09022 
1552  0(2)=1 .42491 

1554  0(3)=. 84331 

1555  0(4)=-1. 21025 
1554  D(5)=. 30744 

1 558  Ii(4)  = -1. 89355 
1540  0(71=1.84787 
1542  0(8)=. 04591 

1544  0(9)=-. 59040 

1545  D( 10)=-. 70444 
1547  0(111=4.29304 
1570  0(12)=-3. 32708 
1575  GOTO  1405 

1580  0(1 )=D(3)=O(4)=0 
1582  0(2)=-. 14447 

1584  0(5)=. 025 

1585  0(4)=-. 08333 
1587  0(7)=. 05247 

1590  O(8)=D(9)=D(10)=0 
1595  0(11 )=D(12)=0 
1405  H2=S0R(M1 ) 

1410  U=D(1)+R1»D(2)  +(R2*D(3)+R1t2*D(4))  +R3*0(5) 

1415  U=U+(R1»R2«D(4)+R1t3»0(7))  +( R4*D( 8) +R2t2»D( 9) ) 

1420  U=U+(R1+R3*O(10)+R1t2*R2*O(1 1 )+RU4«0(  12)) 

1425  T=E(1)+S0R(K(2))*U/H2 
1430  RETURN 


FIGURE  11  SUBROUTINE  FOR  CORNISH  - FISHER  EXPANSION 


